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Abstract. We perfom a thorough computation of the one- loop corrections from both 
scalar and tensor degrees of freedom to the power spectrum of curvature fluctuations for 
non-canonical Lagrangians in single-field inflation. We consider models characterized 
by a small sound speed c^, which produce large non-Gaussianities. As expected, the 
corrections turn out to be inversely proportional to powers of c^; evaluating their 
amplitudes it is then possible to derive some theoretical bounds on the sound speed 
by requesting the conditions necessary for perturbation theory to hold. 



1. Introduction 

The Cosmic Microwave Background (CMB) represents one of the most important sources 
of information for cosmology. CMB fluctuations from the last scattering surface contain 
precious indications about early Universe physics. The Wilkinson Microwave Anisotropy 
Probe (WMAP) and the Planck satellite have been helping collect all these "signals" 
from the ealy Universe and with increasing accuracy and sensitivity [1, 2, 3, 4, 5, 6]. 
Observations have confirmed that many of the existing models share a satisfactory agree- 
ment with experiments. The new data are expected to further help understanding what 
is the correct model. Great expectations rely on non-Gaussianity [7, 8] and on the pos- 
sibility to improve the measurements of the two-point function of cosmic fluctuations. 
This has encouraged the theoretical work on three- and four-point correlation functions 
of the cosmological fluctuations and on the study of loop corrections to the tree-level 
correlators as well. Both the former and the latter are in fact deeply mo del- dependent 
since they are related to the specific interactions characterizing the primordial fields. 

Loop corrections have received a renewed interest since the works [9, 10], where one 
of the main issues to be addressed is whether loop corrections can carry information 
about the history of inflation and whether they can become large. Among the various 
results is that the power spectrum of the curvature perturbations acquires a new (com- 
pared to the tree-level) dependence from the external momentum, i.e. of a logarithmic 
kind (~ Ink), arising from ultraviolet divergences (see also [11, 12, 13, 14, 15, 16, 17]). In 
fact, recently the authors of Ref. [18] showed that such a dependence must be computed 
with care, considering a cutoff of the order a{tk)A, where is the time of horizon exit for 
the mode k during inflation, so that the logarithmic dependence ln{k/a(tk)A) turns out 
to be actually fixed by ln{H / h) |. However, loop corrections to the cosmological corre- 
lators have attracted a lot of interest mostly because of the following reason: they are 
characterized by infrared divergences (~ InkL, where is an infrared cutoff), and by 
the presence of secular terms growing like the logarithm of the scale factor (~ ln{—kri)), 
which diverge in the limit of late times. A debate has then been stimulated about the 
significance of these divergences. 

As to the InkL terms, an interpretation for them has been proposed in terms of the so 
called "small" and "large" boxes: the former corresponds to considering an (inverse) 
infrared cutoff L of the size of our observable Universe, which also implies that k is not 
much smaller than L~^; on the other hand, a "large" box corresponds to choosing a 
cutoff L ^ k~^ . When dealing with observable quantities, the natural cutoff should be 
the "small" box, whereas any correlator computed considering the "large" is expected 
to be interpreted as some quantity that is unaccessible by measurements [19, 20] (see 
also [21]). In this view, the quantities ^ InkL are not large and therefore should not 

X Notice that the resuhs of [11, 12, 13, 14, 15, 16, 17] have been obtained using different regularization 
techniques. In particular, in [15] we employed a cutoff on comoving momenta, rather than resorting to 
dimensional regularization. 
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be interpreted as actual divergences. Another approach might consist in a resummation 
of these divergences through Renormahzation Group techniques (see the interesting ex- 
ample of Ref. [22]). 

As to the ln{—krf) terms, their presence was also discussed in [9, 10], where it was 
pointed out that these terms are generally expected in the loop corrections, unlike any 
positive power law dependence from the scale factor. This divergence appears to be 
fictional since observations are "made at a finite time" [14], i.e. rj ^ 0. Regarding this 
issue, the SN formalism [23, 24, 25] offers a nice way out because two different times are 
involved: an initial one, at which the correlators of the scalar field is evaluated, an a 
final one, which corresponds to the time of observations. The final time is a late time, 
the initial can be arbitrarily chosen and a very convenient choice is to set it just a few 
e-foldings after horizon crossing. This means that the quantity \kri\ is expected to be 
quite close to one. 

When it comes to computing resummations of loop corrections, an interesting discussion 
is provided by [27], where a "Dynamical Renormahzation Group" (DRG) technique is 
developed in order to support the perturbative approach at late times (this technique is 
also employed in order to improve the IR behaviour by resumming late time correlators) . 
See also [28] for a recent review on inflationary loop correction. 

Loop correction have been computed for the power spectrum of curvature fluctua- 
tions in the basic single- field slow- roll model: in [13, 14], these corrections were computed 
considering the fluctuations of the background metric, setting however the tensor modes 
to zero for simplicity; in [15] the tensor modes were fully accounted for and a complete 
computation was carried out in order to evaluate the nature of the tensor mode loop cor- 
rections, which were proved to be characterized by infrared logarithmic divergences and 
to be of the same order of magnitude as the scalar corrections. This showed that a full 
computation of one-loop corrections in standard single-field slow-roll inflation requires 
the knowledge of the tensor contributions §. Loop corrections have also been com- 
puted in more general theories of inflation such as canonical multifield [45] and multiple 
Dirac-Born-Infeld models (DBI) [29] (where a multifield DBI model was considered and 
the loop corrections arising from cross-correlations between adiabatic and non-adiabatic 
modes were computed for the scalar fields power spectrum). The latter models are par- 
ticularly interesting because they are part of a larger class of models characterized by 
an inflaton lagrangian that is a generic function P{X, (p) of the inflaton field and its 
first derivatives, X = —^g'^'^dfj_(f)di,(f), with a sound speed Cg which can be small. In these 
models the interaction Hamiltonian is generally characterized by terms that are propor- 
tional to inverse powers of Cs and that can be therefore responsible for non-negligible 
amplitude of the three- and four-point functions [30, 31, 32, 33, 34, 35, 36, 37], as well as 
loop corrections to the two-point function. For models of inflation that comprise higher 
derivatives of the inflaton field in the Lagrangian see also [38, 39, 6, 40]. For more works 

§ Notice that in Eq. (62) of [15] the numerical coefficient muhiplying the In kL term has been corrected 
in [28]. In particular /2 is equal to a factor 2. We thank David Secry for pointing this out. 



3 



on inflationary loop corrections see also [41, 42, 43, 44, 45, 46, 18]. 

In this paper we are going to study single-field P{X, 0) models of inflation and 
formulate a complete computation of the loop corrections to the power spectrum of the 
curvature perturbations, both in the spirit of [15], i.e. by considering all of the met- 
ric fluctuations (including tensor modes), and in the spirit of [47], i.e. using the loop 
computation and the perturbative expansion conditions in an attempt to extract some 
constraints for the parameters of the theory and for the sound speed in particular. 

The paper is organized as follows: in Sec. 2 we list the background equations 
of motion and derive some fundamental expressions for the parameters of the theory 
(summarized in Table 1); in Sec. 3 we recall the 6N formalism, with some comments on 
its applicability to generic P{X, 0) inflation models, and we employ it for the calculation 
of the power spectrum of (; in Sec. 4 we review the main steps for deriving the interaction 
Hamiltonian and report the leading order terms that will enter the loop calculations; in 
Sec. 5 we report the order of magnitude of the scalar and of the tensor loop diagrams; 
in Sec. 6 we present our flnal results; flnally, in Sec. 7 we draw our conclusions. The 
calculations we performed are quite lengthy, so many of the details are left to the 
Appendices: in Appendix A we solve the Hamiltonian and momentum constraints for 
P{X, (p) models; in Appendix B we present the complete expressions for the third and 
fourth order Lagrangians, including tensor modes; in Appendices C and D we review 
the computation of, respectively, two-vertex and one-vertex loop diagrams; flnally in 
Appendix E we provide some of the polarization tensor equations that are needed for 
tensor loop calculations. 



2. Background equations of motion 

We will consider a class of theories described the following action 

S=^J d^x^^ [mlR + 2P(X, 0)] (1) 

where R is the four dimensional Ricci scalar, mp = (SttG)"^/^ and P is a generic function 
of the inflaton fleld and its flrst derivatives X = —}^g^'-^d^(f)dy(j). 
The background equations are [31, 49] 

2H + = -P, (2) 
H' = ^ [2XPx - P] , (3) 

X {Px + 2XPxx) + 2V3{2XPx - Pf'^XPx = V2X (P^ - 2XPx^) , (4) 

where a dot indicates a derivative w.r.t. cosmic time, P^ = dP/d(f) and we deflne 
X = 0^/2 (we have set mp = 1). In the case of standard single-fleld models of slow- roll 
inflation P = X — V{(j)) where V^(0) is the potential for the scalar fleld. 



4 



Let us introduce some "slow-variation" parameters 

which are assumed to be smaller than one, as for the slow-roll parameters in the standard 
case. However, in P{X,(j)) theories it is not generally correct to talk about "slow-roll"; 
in fact, P is only known to be a generic function of X and (p, therefore the smallness of 
e and r] does not necessarily indicate that 0^ <^ and -C \3H(f)\. 
The speed of sound in these models has the following expression 

Px + 2XPxx- ^ ^ 

We do not assign a specific value to Cg at this stage, the only assumption we make is 
that its derivative is small, more precisely 



which appears as a natural choice dictated by the expression of the spectral index of 
scalar perturbations for these models (n^ — 1 = — 2e — r] — s, see, e.g., Ref. [50, 31]). 
The three parameters e, rj and s in P{X, 0) theories are dubbed as "slow variation" 
parameters because e and t] indicate that the variation of H w.r.t. time is quite slow, 
i.e. the expansion is quasi de-Sitter; similarly, a small s indicates a slowly varying sound 
speed. 



2.1. Useful relations in terms of slow-variation parameters 

A computation of the inflaton power spectrum up to one-loop requires to expand the 
Lagrangian up to fourth-order in the perturbations. The perturbative expansion of the 
Lagrangian is expected to be structured as the sum of a very large number of terms, 
each one weighted by coefficients given by the zeroth order part of P and its derivatives 
w.r.t X and 0. It is then useful to find out the dependence of these coefficients from 
the "slow- variation" parameters and from the sound speed Cg. The main results of this 
section are summarized in Table 1 where we give the order in slow-variation parameters 
(and in the sound speed) of the various coefficients. 

To begin with, the expression of XPx easily follows from the combination of (2) and 
(3) 

XPx = -H, (8) 

which can be plugged in (5) to have 

XPx = eH"^ (9) 

(as a check, in the standard case Px = 1 and the equation above gives e = (0^) / (2H'^)). 
The expression of X^Pxx in terms of slow variation parameters and of the sound speed 
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is derived just as easily from the combination of Eqs. (6) and (9), i.e. 

X'Pxx = I - l) (10) 

The expressions for the higher derivatives of P can be obtained by successive derivations 
of (9) and (10) w.r.t X and 0. We thus expect a dependence not only from the 
parameters e, r], s and Cg, but also from their derivatives. 

We will now derive some equations relating the derivatives of P with each other and with 
Cg and the slow- variation parameters. It is useful to define the following combinations 
of derivatives [31] 

J: = XPx + 2X^Pxx, (11) 
A = X'Pxx + ^X'Pxxx, (12) 

n = X Pxxx + tX Pxxxx- (13) 
5 

From (9) and (10), we get 



e 



H\ (14) 



There is no such simple expression for A and 11 in the general case. The former, for 
instance, can be written as 



A 



where we introduced the partial derivatives of the sound speed w.r.t to X and 
as Cs = Xcsx + (pCsrf). In order to derive Eq.(15), one has to first take the derivative of 
Eq. (6) w.r.t. X, then multiply both sides by X and use the definition of A in (12) and 
of the sound speed. 

It will also be convenient for future calculations to decompose e as the sum (here 
Ex is not the partial derivative of e w.r.t. X but merely the definition provided in (17) 

e = e^ + ex, (16) 

where 

^^^~lP^' ^""^'IPdX' ^ ' 

In standard single-field models of inflation E, A and 11 are equal to zero, whereas and 
Ex become 

where = dV/dcj) and 

1 • •■ 1 00 
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Summing up and Sx from (18) and (19), we recover e = (0^/2)(l/if^). If the slow 
roll condition <^ \SH(p\ is applied, ex is evidently subdominant w.r.t. e^. 
For a general P{X, 0) theory one finds 



3(1 + c2) 
Equivalently we have || 



3(1 + C2 - AXCsCsx] 



(jyXPx^ 



2[e + s) — Tj ^ TTTT-^ — 2' 



He. 



Eq. (20) can be derived by using Eqs. (3), (4), (17) and noticing that (3) implies 

dH 

- 2XPx4, = -QH—, 

so that 

e 

The previous equation can also be written as 
2XHc't 



by using X = cfxp and 
■; d 



r]-2s 



HP 



X 



dt Px 



2e 



r]H 



sH 



+ XPxx 



2P 



X 



(20) 

(21) 

(22) 
(23) 

(24) 
(25) 



Eq. (20) follows from equating (23) and (24). In order to derive Eq. (21), we simply 
take the time derivative of S in Eq. (11), i.e. 

t = XEx + (26) 

and we use the definition of the sound speed, together with Eq. (23). It is easy to verify 

that (20) and (21) agree with each other. 

Other useful and easily verifiable equations are the following 



2i/3 

^XPx^ 
H^ 

^X^Pxx^ _ e 



35 - 2e'^ + 61] + ?>excl 



eT]-2e^ + ?,ex{l + c]) 



(27) 
(28) 



H^ 



, ^^x , 2 i\ , 9A ^ex 



2 e 



\{cl-l) (r/-2e + 3^(l + c; 



2 



(29) 



II E.D. thanks Sarah Shandera for corresponding about the contents of this section and for spotting 
an error in the first version of these calculations. 
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2.1.1. Some worked examples 

In order to get a hint of the typical values for these quantities in some of the known 
models, we will derive the values of Sx and Csx for three cases: 1) DBI model; 2) 
P{X,(j)) = f{X) + g{(j)) models (where / and g are generic functions of their ar- 
guments), of which the canonical slow- roll Lagrangian represents a special case; 3) 
P{X, (p) = f{X)g{(()) models (of which K-inflation is a special case). 

• DBI model. 



P fi^ v/l-2X/(0) , 1 

PnniiX, 0) = + - r(0). (30) 

In this model we have 

cl = l- 2Xf (31) 



and 



so 



and 



Px = -, (32) 

Cs 

(33) 



1 - c2 1 - 

^ = ~2X^ ^ 2i^- ^^^^ 
Let us compute the order of ex- From Eq. (20), we know that this can be derived 
by computing 0XPx</.- Taking the derivative of Px in Eq. (32) w.r.t. the scalar 
field we get 

^^P^^ = (1 - 2'x/)3/2 = ^ • ^35) 



Replacing the derivative of Eq. (34) w.r.t. time 
Eq. (35) becomes 



<l>XPx^ = r—^ ^ -^H' (37) 

L.s Cg 



where, in the last line, we assumed that rj ~ 0{e) ~ 0{s) and that < 1 (we 
are interested in small sound speed values). Replacing this result in (20) we get 
Ex — 0{e'^/c^). As a double-check, we can calculate the expression of ex by taking 
the time derivative of (33) and using the definition of 77; the result is 

= {v + s-2e). (38) 
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Let us now evaluate Cgx- Taking the derivative of (31) w.r.t. X we get 

2csCsx = -2/ = , (39) 

then XcsCsx — 0{l) if < 1. 

P{X,ct>) = f{X)+g{cP) models. 

In these models Px4> = 0, therefore ex — Oi^e^) from Eq. (20). The sound speed is 
independent of 

JX + ^ A J XX 

so from = Xcsx and using Eq. (20) and (23) we get 

2Xc.c.x = (l + c^) . 
This brings Xc,c,x ^ 0{l) if s ~ 0(e). 

P(X,0) = /(X)(7(0) models. 

The order of Sx can be computed knowing that 

0XPx0 = eH^^-, (42) 

so 

^x = ;77T^f25-r/ + ^V (43) 



(41) 



3(1 + c2) V ' gH^ 

The first of Eqs. (41) also applies for these models, so XcgCgx — ^s/sx- In a 
subclass of K-inflation models, known as "power-law K-inflation" [51, 50] 



where 7 is a constant. The sound speed in these models is given by 

„2 7 



(44) 



(45) 



' 8-37' 

so the regime of small corresponds to considering small values of 7. In this regime, 
it is possible to show that g/ gH ^ 7, from which we have Sx — exO{e+c^g) ~ 0{e^) 
(it is easy to verify that e ^ ~ 7 for (44)). 

We now know what are the typical values of Sx and Csx in some of the known P{X, 0) 
models. Let us see how this information can be used to derive the orders of magnitude 
of the remaining derivatives of P. 

Let us begin with the derivatives w.r.t. X. From the definition of the sound speed we 
have 

A-P,, = £^ (^ _ 1) , (46) 
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so 



This can be equated to 
d 



[XPx + X^Pxx) + 



XP 



X 



X^P 



XXX 



In a similar way, we can derive 



X^P 



xxxx 



— 3X Pxxx + 



+ 



1 \ 2XCsCsX 2X'^CsCsXX 



-2XCsCsX 



^dX ^^'^'^^^^ ~ + X^Pxxx, 

in order to solve for Pxxx- Using e = XPx/H'^, we get 



2XCsCsX 



(47) 



(48) 



(49) 



2XCsCsX 



(50) 



where Csxx = d'^Cg/dX'^. Notice that the equations above are valid for a generic P{X, (p). 
The remaining derivatives can be obtained in a similar way. For instance, we know from 
Eq. (27) the expression of (pP^ as a function of e, t] and Sx- By taking the derivative 
of (27) w.r.t. and using (28), it is possible to easily derive the order of magnitude 
of in terms of the flow parameters e, 77 and s, ex and the first derivatives of 77 and Sx- 



For practical reasons, given the large number of terms in the Lagrangian, we will fix the 
order of magnitude in terms of e, s, Sx and Csx- In the case where 

Sx^O f 3^ , (51) 



and 



XCgCsX ^ 0(1), 



0{e) 



(52) 



we include all DBI models, see first column of Table 1. 
On the other hand, when 



ex ^ 0{e^) (53) 

together with the condition in Eq. (52) we manage to include all P{X, 0) = f{X) +g{(f)) 
class of models (notice that Eq. (53) implies ix — 0{e^), ex — 0{e'^)). We report the 
results for all these cases in the second column of Table 1. 

Among the P(X, 0) = f{X) x g[(f)) models, we list the orders of magnitude for the 
case described by Eq. (44) in the third column of the table. Notice that, as shown in 
this section, the coefficients in Table 1 are not all independent from one another and 
can be combined in a smaller number of independent parameters, such as in Eqs. (11) 
through (13) (see e.g. [31]). Another way of defining these parameters is presented 
in [39, 6] whithin an effective field theory framework, where the ones that are linearly 
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Table 1. Order of magnitude of the partial derivatives of P for three classes of 
P{X,(t)) models. Class I is identified by Eqs. (51) and (52) and includes DBI-like 
models. Class II is identified by Eqs. (53) and (52). Class III corresponds to the 
model described by Eq. (44). 





class I 


class II 


TTT 

class ill 


XPx 




/ ^\ zj2 


/' \ zj2 


V2 D 


{e/ci)H' 


1 ^ 1 2 \ TJ2 


zj2 


J^XXX 


/ / 4\ tt2 


/ / 4\ tt2 





^ ^XXXX 


f / 6\ tt2 


{e/c°)H- 





J. J ) 


[e)H 


[e )H 


















(^3)^5 


4>'^ P(i)<j,(i)<j, 















(7^)^' 


<pX^Pxx^ 









4>X^PxxX(f> 










(j)'^X Pxtjxj, 







(7^)^' 


(j)^X Pxt/xj),/) 









(j)'^X'^Pxx<t>4> 










independent are reported to leading order in slow-roll. 

In our loop calculations, we will consider the leading order terms in the Lagrangian, 
restricting to the cases included in the first two columns of Table 1. Notice that, tak- 
ing the leading order terms only, i.e. neglecting the coefficients with derivatives w.r.t. 
(j), allows to consider at the same time both class I and class II models on the same level. 



3. (5N formalism applied to P{X, (j)) models 

The (5N formalism [23, 24, 25, 26] relates the comoving curvature fluctuations defined 
on a uniform density temporal slice at time t to the initial perturbations of all the fields 
defined on a fiat slice at time t* < t. In standard slow-roll inflation, we have 

„ /V(")(t t ) 

as.t) = Y, , mx)^)\ (54) 

n 

where, N^'^\t^t^) = d^N/dip^ represents the n-th derivative of the unperturbed number 
of e-foldings N{t,t^) w.r.t. the unperturbed values of the scalar field at time t*. The 
computation of the power spectrum of the curvature perturbation produced during 
inflation can thus be made in two steps. First one needs to compute the power spectrum 
of the initial inflaton field perturbations 50*, and then the 5N formalism (54) can be 
employed to obtain the final value of the power spcetrum at a later time t. As discussed 
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in Refs. [25] and [15], even in the case where the tensor fluctuations are taken into 
account, they will not (explicitly) appear in the previous expansion. In the case of 
general single-field models of inflation it is convenient to choose the initial time a 
few e-folds after the crossing of the "acoustic horizon" which goes as Cgf]- A priori not 
only the field S(j) but also its first time derivatives would in principle appear in the 5N 
expansion, the reason being that the field (j) obeys a second order differential equation 
and therefore its background trajectory can be uniquely determined once the field and 
its first derivatives are specified on the initial slice. In the case of standard single-field 
models, assuming the usual slow-roll approximation, such as 3H(j) ~ —V'((f)), implies 
that (f) and <p are no more independent from each other, therefore the initial state for 
the system is completely determined by (p and it is not necessary to also include 6(f) in 



Since for the more general class of models with non-canonical kinetic terms in the 
Lagrangian P{X, 0) inflation can be attained even when the usual slow-roll parameters 
are not small, one could wonder whether such a simplification can be still employed, so 
as to still make use of the simple expression (54). For P{X,(j)) models, the equations 
of motion are in fact more complex. Nevertheless, the smallness of the slow-variation 
parameters and of allow to get rid of in (4). Both for P = f{X) + g{(f)) and for 
DBI models, we have indeed verified that (4) reduces to 



by taking into account the order of magnitude of the different terms appearing in Eq. (4) 
and the results of Sec. 2 (see also [48] for another derivation of Eq. (55)). Like in the 
standard case, this is an equation relating and its first time derivative, which are, 
again, no more independent (it is easy to check that a similar result can be obtained 
for "power-law" K- inflation as in (44)). In Ref. [34] there is a hint about the use of 
the 6N formula for general single-field models of inflation and the reader can deduce an 
alternative way which validates the use of Eq. (54). 

We will be using Eq. (54) in the computation of the power spectrum P^ of the curvature 
fluctuations 



(54). 




(55) 



(56) 



Up to one-loop it can be expressed as [44, 14] 




(57) 
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In Eq. (57) Ptrcc(^) is the tree level power spectrum of the inflaton field [50] 

{^<f>^MkJ* = {2nfS^'Hh + k^2)^^^, , (58) 

with the Hubble parameter, and also the other quantities, evaluated at horizon exit 
(when kcs = aH). -PonG-ioop(^) is the one-loop contribution to the powers spectrum of 
the initial inflaton fluctuations 

-Pone— loop (^) -Pscalar(^) ~l~ -Ptensor(^) ; (59) 

where Pscaiar, accounts for the contributions coming from scalar loops, whereas -Ptensor(fc) 
is due to loop diagrams where tensor modes are involved. Both of these corrections will 
be computed in this paper. 

Finally, the second line of (57) includes the integral of B^{ki, k2, k^), the bispectrum 
of the scalar field (see [31] for the evaluation of for general P{X,(j)) models). Since 
(Ar(i))" = {Mp^){Px/2e) the tree-level power spectrum of the curvature perturbation 
reads 

where all the quantities on the r.h.s. are meant to be evaluated at horizon-crossing. 



4. Perturbative expansion of the action 



The power spectrum of 6(f) can be computed using the Schwinger-Keldysh (also dubbed 
as "in-in") formula 

r pv 

(fi|50,;;(r/)50,;;(r/)|fi)i^ = z(oT 50,-;(r/)50,-;(r/) / rfr/' fif+,(V) - ^rn*(V)) 0) (61) 



+ 



T 



J — oo ' . 

2 (O|T[50,-;(r/)50,-;(r/) r dri (^HUr]') - Hr,{7^' 



dv" (H^ntiv") - Hi^tiv" 



where Hint is the interaction Hamiltonian and the plus and minus signs indicate the 
propagators in the in-in formalism [52, 53, 54] (see also [9] for a complete review). 
The field operators appearing on the right-hand side of Eq.(61) can be expanded in 
terms of their Fourier modes 



6(f){x,r]) = / d^ke'"'' \aj:6(f)k{v) + a^Mliv) 



with creation and annihilation operators given by 

{27rf5^^\k-k'), 

= {27r)'6^^\k-k')6,.. 
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The eigenf unctions for scalar modes are ^ 

6Mv) = /tAtts (1 + ^^^^^) (62) 
whereas for the tensor modes 

7.(r/) = -^(l + zMe-^'''. (63) 

One loop corrections to the power spectrum are sourced both by three and four-leg 
interactions, which means that the perturbative expansion of Hint should be carried out 
up to fourth order in the field fluctuations. 

It is convenient to adopt the 3 + 1 Arnowitt-Deser-Misner (ADM) splitting for the metric 
tensor. In the spatially fiat gauge the metric is 

ds"^ = -N'^dt^ + h,j {dx' + N'dt) {dx^ + Wdt) , (64) 

hij = a\t){e^)ij, (65) 

where a{t) is the unpertrubed scale factor, is a tensor perturbation with di'~^ij = 
'jii = (traceless and divergenceless) and det(e'>')ij = 1. and A^j are called "lapse" 
and "shift" functions and can be seen as Lagrange multipliers: they obey first order 
differential equations (respectively called the "Hamiltonian" and the "momentum" 
constraints) that can be solved to derive the expression of A^ and Ni as functions of 
the field perturbations 6(j) and 7^. The action (1) can be written as follows [55] 



S=^J dtd^xVh [NR^ + 2NP + N"^ {EijE'^ - E^)] (66) 
E^J = I {h^, - V^NJ - VjN,) , (67) 



where 

E.. = 

2 

E = h'^Eij, (68) 

R^^^ = -^daabdaab. (69) 

Following the general procedure of Refs. [55, 31], the next step is to write and solve 
the constraint equations, replace them in the action and expand the Lagrangian up to 
fourth order in the perturbations 6(j) and 7jj. We solved the constraints in [15] (see also 
Appendix A of this paper for a review of the computation). The complete expressions 
for the third and fourth order action instead represent a new result and are provided in 
Appendix B. 



% Notice that, at least for the cases we are taking into account, the prefactor l/Vcs-Px in Eq. (62) 
has a time derivative that is much smaller compared to the quantity itself (and even exactly equal to 
zero for the DBI models discussed in Sec. 2.1.1); this is important in order to be allowed to treat this 
prefactor as a constant when integrating over time. The same applies to the (f/>)~" {n > 0) factors that 
appear in the interaction terms of the action (see for instance Eqs. (70)-(71)). 

A perturbative expansion of the Lagrangians for P{X, 0) models was done in [34] , but with a different 
approach as far as the tensor modes are concerned; here we expand the tensor mode fluctuations in a 
similar way as we do for the scalar ones. 
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4.1. Leading terms in the third and fourth order actions 



Before we proceed with selecting the leading order terms according to the assumptions 
we made in Sec. 2 (i.e. according to columns I and II of Table 1 therein), it is important 
to make a preliminary remark about the diagrams with tensor loops to fourth order. 
As shown in [15] and also pointed out in [45], the momentum integrand functions for 
these diagrams are expected to be independent of the external momentum k, therefore 
they can only exhibit renormalizable power-law ultraviolet divergences. Their final 
contribution should then be included in a leftover renormalization constant. For this 
reason, the computation of these contributions will not be performed. We have followed 
the same procedure to account for the one-vertex scalar loop diagrams (see below and 
Appendix D for a more detailed discussion of this point). 

With this remark in mind, we will report in this section the action to third and fourth 
order. For the third order action, we report the leading order contributions from 
the scalar interactions and all of the contributions from the terms involving tensors 
(determining the order of magnitude of the tensor contributions is not as trivial as for 
the scalar ones, so we decide to report all of them at this point) 

S • P / ■ \ 

+ ——jiijdklijdkd~'^5(l) + y - (I^Hdiljkdiijk + 8H%jdi5(j)dj6(j)j 



-^j^iij6(f)didjd~''^6<i) - j^dajkdkd~'^6(j)didjd~'^c 



(70) 



For the fourth order action, we report the leading order scalar terms only 
J 2 [ 12 ' (f)^H (p^H V ' / V ' 

— a2^4> H — 7^ (diScj)) a2 + 6X6(t)di6(f)di92 



iH2 



(71) 



We indicate [diScj))^ = SijdiS(j)dj6(j) and sums are intended over repeated spatial indices. 
Also, we define 



02 = d^' (72) 

d2 = d-^ [A + B] (73) 

= [d^^dj (A + 5) - C] (74) 
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where A = {d'^S(f)) 50, B = Sij {diS(f)) (djS^j and C = (6»/0) 50. 

The complete (i.e. without any slow-variation parameter approximation and including 
all tensor mode interactions also at fourth order) expression of the action up to fourth 
order were computed and can be found in Appendix B. It is important to point out that 
the selection criteria we adopted in order to identify the leading scalar contributions, al- 
though they apply to a large number of P{X, 0) models (as described in Sec. 2), are not 
valid for all of them. However, our complete expression of the interaction Lagrangian 
offers the possibility of computing loop corrections (as well as higher order correlators) 
in P{X, 0) models other than the ones we considered. 



4-.2. Calculation of the interaction hamiltonian 

In a theory, like the one we are dealing with, where time derivatives of the fields appear 
in the non-quadratic terms of the Lagrangian, the relation Hint = —Li between the 
interaction Hamiltonian and the Lagrangian does not generally hold [56, 13, 32] (also 
see [15] for a quick discussion on this). 

We follow the prescription reviewed in [45] for computing Hint- It is convenient to 
schematically write the Lagrangian up to fourth order as follows 

Gi"^50' + •"•''^^ + jf''^''^^ + C^<50' + + /if ) 502 + /if ^50 + jM , (75) 



L 



where /, g, G, j and h are functions of the arguments (and their derivatives) that are 
indicated in their upper indices, whereas the lower index refers to their pertubation 
order. They can be easily read from Eqs. (70) and (71). Notice that we have for 
simplicity ignored the tensor indices. 

The momentum densities for the scalar and tensor perturbations are 

(77) 

from which we have 

50 = 50(1) + 50(2) + 50(3) , (78) 



where 



- ^ (79) 

,(2)_ gj'"' Sgl'^\('^)' G(^)vr(^^) ^(^)vr(^) ^^^^ 



2)(^^) 2(/(^^))' 
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1(3) 



9 Uo 



) 



16 (/o^ 



^(7)^(2) G'S^^50(2) 



4 /o^ 



4 /, 



2 /, 



Jo 



and 



where 



7 



^(1) +^(2) +^{3) 



(1) 



vr 



(7) 



7 



^(3) 



(7) 



(7) 
92 

p(7) 



2/o 

" 2/(-) 



, ^((50) ^(7) 
^W)^(2) 



^W)^(7) 



2 (/o^^^ 



2 ' 



Jo 



(7) 



The Hamiltonian density 



(81) 

(82) 

(83) 
(84) 

(85) 
(86) 



is derived plugging in Eqs. (78) and (82) into its definition. We finally replace the 
conjugate momenta with the time derivatives of the field perturbations as derived from 
the free Hamiltonian density Hq 



7 

and the result is 



(3) 



(87) 
(88) 

(89) 



for the third order part of the interaction Hamiltonian and 



2h. 



h 



+ 



92 



^Po 



7 



AH) , ^ 
J4 + ^ 



4(.f 



9i 



\ 



Jo 



(7) 



^Jo Jo 



Jo 



(7) 



(50) 



AH) 



+ 2/0^^^)^^ 



i^^+.2 



(7) 



+ 



AH), 



(7). 



^2 



{H)\ 



J 



+ f^^ + 9^2 



(7) 



4 /d'^'^'' fo''^ 



(90) 
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for the fourth order part. In the case where the tensor perturbations are set to zero one 
can verify that Eqs. (89) and (90) reproduce Eq.(13) and (14) in [32]. 

5. Calculation of the loop corrections 

In this section we will provide the order of magnitude of the one-loop corrections arising 
from the third and the fourth order interaction Hamiltonian. 

The order of magnitude of the different diagrams can be quickly estimated for the di- 
agrams with scalar loops. In fact, for these diagrams it is possible to factor out of the 
time and of the momentum integrals all of the if, c^, Px, s and all other parameters 
of the theory. For the diagrams with tensor loops it is instead necessary to evaluate 
the time integrals before the overall coefficient can be determined. We first explain 
the prescription for deriving the order of magnitude of the final result for each type of 
diagram and, after that, we provide the actual orders of magnitude for all the leading 
order diagrams. We report their precise expression in Sec. 6. 

• Scalar loop diagrams 

— Two-vertex diagrams 

Let us consider two three-scalar interaction verteces, which we will indicate 
as Vi and V2, from the interaction Hamiltonian (i.e. / d^xdrja'^H^^^ D 
J d^xdria'^{ri){Vi + V2)). The verteces Vi {i = 1,2) will in general be expressed 
as Vi = Vi{H,e, Px, ■■■)fi{S(f>,d^6(f)), i.e. as the product of a function of 
some parameters of the theory, multiplied by the scalar field fiuctuations and 
their derivatives. The prescription for calculating the overall (dimensionless) 
coefficient C of the diagram built from Vi and V2 (as in Fig. 1) is the following 

C ^ (v^) X (v,) X (^7^)' X (-^y X (Hr^^^-, (91) 
the complete analytic expression of the diagram being given by 

(%,(^*)%,(^*)) 3 + h)CG{x*)F{k), (92) 

where F is a function of the external momentum k = \ki\ = \k2\ with 
dimensions (mass)~^ and G is a function of x* = —kCsTj*. The coefficient 
{H/ (y/csPx))^ comes from having eight eigenfunctions like the one in Eq. (62); 
(Cg/if^)^ comes from having two temporal integrations drja'^l't]), which we 
rewrite in terms of the dimensionless variable x as {dx/x^){—k^c^/H'^), using 
a ^ {Ht])"^, valid for a quasi de-Sitter spacetime; finally, extra factors of H 
are needed for each temporal derivative (spatial derivatives do not contribute 
to determine C) appearing in the vertices, as we can see from the following 
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Figure 1. One- loop two- vertex corrections to Pscj, from scalar modes. 



Figure 2. One-loop one-vertex corrections to P^^ from scalar modes. 



Figure 3. One-loop tensor corrections to Pscf, from two-scalar one-graviton 7 {Scj)) ) 
interactions. The continuous line represent a scalar propagator, the dotted line a tensor 
propagator. 



expression 

SMV) = -^f^^^'e-. (93) 



— One- vertex diagrams 

If we consider diagrams with only one vertex V = v{H,e,Px, ■■■)f{^4>,dfj^6(j)) 
arising from scalar interactions (see Fig. 2), the prescription will be slightly 
different 

given that the number of time integrals has now gone from two to one and the 
number of eigenfunctions from eight to six. 



• Tensor loop diagrams 

The evaluation of C for diagrams with tensor modes running in the loop is not 
as straightforward as for the scalar diagrams because the sound speed cannot be 
so easily factored out before actually solving the time integrals. In fact, the time 
integrals will be over a combination of derivatives of the eigenfunctions of the scalar 
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Figure 4. One-loop tensor corrections to Ps^ from two-graviton one-scalar 5(f) (7) ) 
interactions. 



ik{'n) = — 7^—^(^'^'- (95) 



modes from Eq. (93) and of the eigenf unctions of the tensor modes 

2H Hx^ 

However we expect that the one-loop corrections involving the tensor perturbations 
will be subdominat with respect to pure scalar contributions, the reason being that 
the scalar field fluctuations depend on the sound speed, while the tensor ones do 
not, and the order of the loop corrections is dominated by inverse powers of the 
sound speed. Such an expectation is indeed confirmed by the estimates below. 

— Two-scalar one-graviton interactions 

If the vertices Vi and V2 indicate an interaction between two scalar and one 
tensor degrees of freedom (see Fig. 3), the prescription for the C coefficient 
becomes 



C ~ {Vi) X {V2) X (^7==) X (Hf X (^-^^ X ^Hfirneder ^ ^ ^^^^ 

(96) 

where the factor comes from having two tensor eigenfunctions in the 
diagram and (/ dxdx) reminds us that extra factors are expected from the 
time integrals. 

Two-graviton one-scalar interactions 
Finally, for diagrams as in Fig. 4, we have 

C ~ {vi) X {V2) X {^L=\ X {Hf X X (//)»**'-- x( ! rfxda;) . 



(97) 

5.1. Order of magnitude of all leading order diagrams 

We are now ready to list the orders of magnitude of our diagrams. 

Let us begin with the scalar loop diagrams. The leading order three-scalar interaction 

Hamiltonian terms are J d'^ xdrja"^ H^^] D J d^xdrja"^ [Va + H] , where 

Va = (98) 
H = -^7^ ^50, (99) 
i.e. (from Eq. (70)) 

Va = -^Pxxx^'S<p\ (100) 
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xx(, 



(101) 



From the combination of Va and H we get three one-loop diagrams with similar orders 
of magnitude in terms of the parameters of the theory 

j^i 1 



(102) 



where i,j = a,b and, as usual, the Fij are functions of the external momentum with 
dimensions (mass)"'^ (here we reintroduce the Planck mass), while Gij are functions of 
X* = —kcgiri^. The explicit computation of the corresponding one- loop corrections to the 
power spectrum of 5(j) reads 



:,(^m,(^*))3 5^'H^i + ^2)vr 
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[X^Pxxx) 



{X'^Pxx) (X^Pxxx) 



Ml 



Gaa{x*)ln{A/H, 



k^ 



2[X'Pxx)' 



e^clPxMl 



G,,{x*)ln{A/H, 
k^ 



(103) 



where Gaa, Gab and Gbb are functions of x* as provided in Eqs. (C8) through (CIO) and 
A is a fixed physical cutoff that was introduced while integrating over momentum [18]. 



The result in Eq. (103) can be rewritten after expressing Pxx and Pxxx in terms 
of the slow- variation parameters and in terms of the sound speed, according to Table 1. 
From Eq. (46) we have 



X^P 



eH^Mj, ( 1 



where, in the last step, we have as usual assumed that we are in the regime of small 
sound speed. From Eq. (49) and for models which satisfy the constraint XcsCsx ~ C^(l); 
we have 



X^P 



XXX 



eH^Ml 



Using (104) and (105), Eq. (103) becomes 



STT^edPxMf, 



From the fourth order action, the leading order vertices are 

H^^l D J dvd'xa\t) (V(o) + V^i) + V^2) + V(^3) + , 



where 



PxXX^ff^ ~ r -,2 



Vi 



(1) 



24- 



H 
AS 



-a2(- 



3 ; ' 



(105) 

In {A/H,)' 
k^ _ ' 
(106) 

(107) 

(108) 
(109) 
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AE 

= -l2-^dM\ (110) 

V(3) = — ^2^^2 {5ijdi5(j)dj5(p) , (111) 

1^(4) = -Q^^S<P6,, (dM) {9,62) . (112) 

Notice that the leading order vertices in H^^l from Eq. (90) (e.g. ~ h^^'^^S^'^) have not 
been taken into account since they provide disconnected-hke contributions, similarly to 
what happens for the one-vertex diagrams with graviton loops (see also Appendix D for 
an accurate discussion on one- vertex diagrams). 
Using the abbreviations of Sec. 4.2, these terms can be condensed as 

i^S = -hf'^S^' - /^f ^50. (113) 

The order of magnitude of the one-loop diagrams built from the vertices (108) through 
(112) is given by 

(%.(^*)%.(^*))^ = + k2)^-±-^G,{x*)F.{k) (114) 

where now z = 0, 1, ...,4. These are evidently subleading compared to (102). 

Let us now move to the diagrams with tensor loops. The interaction terms with two 
gravitons and one scalar are 

^91 = -::T7Tiijdklijdkd~~^6^, (115) 
2H(p 



= ^^ir^s<f>dajkdajk, (ii6) 

oil CL 



Vg3 = -^iijiij^^)- (117) 

The orders of magnitude of the one loop diagrams resulting from their combinations are 
as follows 

:%,(^*)%.(^*)w = ^^'n^i + ^2 



* 


e 


Ml 


Pxcl 


Ht 

* 


e 


Ml 


Pxct 


Ht 


e 


Ml 


Pxct 


Ht 

* 


e 


Ml 


Pxc's 


Ht 


e 


Ml 


Pxc's 


Ht 

* 


e 


Ml 


Px^s 
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(118) 
(119) 



:%,(^*)%.(^*)).i.3 = 5^'\ki + h)^-^^G,,,,{x*)F,,,.,{k), (120) 

:%.(^*)%.(^*)).... = ^^'H^i + 4) 

:S<PkSv*mM))a.a. = + h)^-^^G,,,,{x*)F,,,,{k) (122) 

IVl p ixCg 

hSv*)S^kM)),s,. = + h)^-^G,,,,{x*)F,,,,{k), (123) 



which were computed after integrating in time. 



Finally, the interaction terms with two scalars and one graviton are 

Vs, = {-i^,di54>dj64>) , (124) 
Vs, = ^ (^iij5<pdid,d~'5^) , (125) 
= ^ (^^^,^~'6^) (da^k) d,d-'6<j>. (126) 

From the combination of these last vertices, we have the following one-loop corrections 

IVl p ixCg 
IVl p JrxCg 

{S<Pj:^{v*)Sh,{vl)s.s, = S^'\h + h)^-^,Gg,s-A^*)Fs,g,{k), (129) 

IVl p rxCg 

(%^(ry*)%^(r/*)),,,, = 5(3)(^i + 4)^^a...(a:*)F,,,,(fc), (130) 
(%,(^*)%^,(^*))...3 = + k2)^-f-,Gs,s,{x*)Fg,g,{k), (131) 

IVl p ixCg 

(%,(^*)%,(^*)).3.3 = + k2)^-^Gg,g,ix*)Fg,g,ik). (132) 

IVl p J^xCg 



As we can see from the previous equations, tensor loop corrections are subdominant 
compared to the scalar ones. The details of the computation of the loop diagrams can 
be found in Appendices C and D. 



6. Final results 



From the expressions of the previous section the leading final result for the one-loop 
correction to the power spectrum of 6(j) turns out to be 

%iv*)S(l)^k,iv*))uoop = (27r)35(3)(fc, + k,) ^ „ ^! ( 



^Ti'^eclPxMl 



9 



7:Gab 
D 



1 ^ In {A/H, 



(133) 



We are now ready to switch from the 5(j) to the C power spectrum. The first line of 
Eq. (57) (i.e. the one-loop "quantum" correction) becomes 



(134) 



{2Txf5^^\h + k2)V^ 



27r2 
k^ 



1 + 



2P^ /8 



-Gaa — -Gab + -Cfofo lu {A/ 
9 O 2 
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where Eq. (58) and the equation (iV^^)) = {Mp'^){Px /2e) were employed (P^ = 
H'^ / {Sir'^ MieCs) , reintroducing the Planck mass). 

The final result for to one loop is obtained by summing up the contribution we just 
computed to the "classical" loops contributions indicated in the second, third and fourth 
lines of Eq. (57), i.e. 

(QWC-^W) D{27rr6^'\h + h)\N('^N('^ I ^i?^(fci,g,|fc,-gl) 



+ \ {N^'^Tj ^,PtrUQ)PtU\h - g1) 

+ iV(l)iV(3)Ptree(A;) / ^Ptroc{q)] ■ (135) 



(27r) = 



The derivatives of are given by 



— 7\/r A/ , ' — i\j2 1 — T\j?, ' l-^'-'^J 



Mp\l e ' Ml' Ml 

where we used the definition of the number of e-foldings N = J Hdt together with the 
results of Sec. 2. Considering Eq. (58) and computing the order of magnitude of P^ *, 
Eq. (135) can be put in the form 

(Cu^MuS^)) ^ {2^)'S^'\h + J) InkL. (137) 

In deriving Eq. (137) the integrals in (135) were solved following the method proposed in 
[57] (the result in (137) is exact up to an 0{1) numerical coefficient); is an infrared 
cutoff that is generally chosen to be comparable to the size of the present cosmological 
horizon (see Sec. 4 of [14] for a complete discussion and for more references on this). 
By looking at the power spectrum of ( with quantum (134) and classical (137) one-loop 
corrections, we can conclude that, logarithmic factors apart, quantum corrections are 
larger than the classical ones. 

Using Eq. (134), it is possible to derive some theoretical bounds on the speed of 
sound, by requesting the one-loop correction not to overcome the tree-level contribution 
for the standard perturbative approach to hold 



2 



where we have neglected a logarithmic term together with other 0{1) numerical 
coefficients from the Gij factors, jl Oii the other hand CMB observations allow a 

* We compute the order of magnitude of tree level diagrams arising from Va and Vb which have been 
found to be the leading order three scalar vertices. 

jj Using Eqs. (C8) through (CIO), the sum appearing in Eq. (134) turns out out to be 
iSGaa/9 - SGab/6 + Gbb/2) ~ 1.16 for ~ 1. 
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measurement of the amplitude of the primordial power spectrum, giving ^ 2.4 x 10 ^, 
so that the corresponding bound for the sound speed is given by 

Cs > 0.9 X 10-2 . (139) 

It is interesting to note that such a bound is very close to the observational constraint 
of Ref. [6], Cs > 1.1 X 10"^ (95% C.L.), derived through a detailed analysis of the CMB 
bispectrum using WMAP 5-year data, ff The more and more refined techniques of 
analysis for the CMB bispectrum used at present and the high precision data available, 
clearly show that a theoretical prediction at the same level of precision is indeed required. 

7. Conclusions 

In this paper we have performed a thorough computation of the power spectrum of 
curvature perturbations ( arising in general single field models of inflation with a non- 
standard kinetic term. In the spirit of Ref. [47], the main motivation comes from the 
fact that in these models large non-Gaussianities can be produced when the sound speed 
Cs is small, with an amplitude of the three-point correlation function (CCC) given by the 
non-linearity parameter /nl ~ c^^ [31]. This suggests that, when computing the two 
point function {(() up to one-loop, strong corrections could arise as well, and invalidate 
the usual perturbation theory approach. A rigorous computation is required to see the 
precise conditions for which the one-loop corrections are under control. In doing so we 
have obtained various results. We have provided the general action at fourth-order (see 
Appendix B), and we have included the tensor fluctuations as well, showing explicitly 
that the loop-corrections involving tensor mode fluctuations are negligible with respect 
to the pure scalar contributions. In deriving the leading order corrections at one-loop to 
the power spectrum of the curvature perturbation, we have expressed the magnitude of 
the various interaction terms appearing in the Lagrangian in terms of slow-roll variation 
parameters and the sound speed, as summarized in Table I. We have applied these 
results to various cases, which include also DBI models. For the usual perturbative 
approach to remain valid, the one-loop corrrections must be smaller than the tree-level 
power spectrum and this imposes an upper bound on the sound speed, Cg > C(l) x 10^^. 
Interestingly enough such a bound turns out to be very close to the observational bound 
derived in Ref. [6]. 
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Appendix A. Solutions for the Hamiltonian and momentum constraints 



In the spatially flat gauge and using the ADM formalism, the perturbed metric is 

ds^ = -N^dt"^ + h,j {dx' + N'dt) {dx^ + N^dt) (Al) 

where 



a\t){e' 



(A2) 



with traceless and divergenceless. The lapse and the shift function can be expanded 
as follows 

N = l + a, (A3) 

N,=d,e + P,, (A4) 

where Pj is divergenceless and a and 6 are scalar fluctuations. Before deriving the 
Hamiltonian and the momentum constraints, we also need to perturbatively expand 
P{X, 0). Let us begin with the expansion of X 

1 



X 



.00^2 



+ h'^ - 



+ 2N"'^N'di5(t) 
di5(t)dj5A = Xo + AX^ 



(A5) 



where = h'^^Nj, Xo is the zeroth order part, i.e. Xq = 0^/2 and AX is the 
perturbation to the desired order. Notice that 0(t, x) = 0o(t)+(50(t, x), but for simplicity 
we will suppress the subscript '0' in the background value of the field (the order at which 
is to be taken will be clear based on the particular context). 

The expressions for the perturbations AXj have been completed as follows: the quantity 
Xo has been factorized 



AXi = 2Xo 
AX2 = Xo 
AX3 = Xo 



50 



ai 



—di6(j)di6(j) 



(A6) 
(A7) 



2ai -Z 



(j) (j) (j) 



(A8) 
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I ycf) J J (f) (f) 


<) 

Up to fourth order we have 

P(X, 0) = P + AX + P^50 + ^PxxAX^ + ^P^^502 + Px0AX50 
+ ipxxx AX^ + ^P<^<^<^50=^ + ^Pxx</,AX250 + ^^P^^^AXdct? 
+ ^PxxxxAX^ + ^P0<^^<^50^ + ^Pxxx^AX^ + ipx^0^AX50=^ 

+ ^PxX00AX2502 

= P + Px (AXi + AX2 + AX3 + AX4) + P^64> 

+ ^Pxx (AXi^ + AXs' + 2AX1AX2 + 2AX1AX3) 

+ ^P^050' + Py0 (AXi + AX2 + AX3) 50 

+ liPxxx (AXi=^ + 3AX12AX2) + ^P^^^6<j)^ 

+ ipxx0 (AXi^ + 2AX1AX2) 50 + ^Px00 (AXi + AX2) 50' 



+ ^PxxxxAXi'' + ^P000050^ + — Pxxx^AXi" 

+ ^Px<^<^<^AXi50=^ + ipxx00AXi2502, (AlO) 

where Px = dxP, P4, = d^P and so on for higher order derivatives (we have for 
simphcity suppressed the subscript '0' on the right hand side in P and its derivatives). 
The kinetic part of the action is exphcitly given to fourth order by 

2XP = 2 (1 + ai + 02) P = 2P(^) + 2aiP(3) + 2a2P^'\ (All) 



where 



P(2) = PxAX2 + ipxx AXi' + ip^<^502 + Px<^AXi50, (A12) 
P(3) = PxAXa + PxxAXiAX2 + Px0AX250 

+ ^PxxxAXi=^ + ^P^0^50' + ^Pxx^AXi250 

+ ipx<^<^AXi50' + ipxxx<^AXi3, (A13) 
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+ ^Pxx (AXs^ + 2AXi AX3) + + ^PxxxAX^AX^ 

+ PxxAX.AX^Scp + ipx<^<^AX2<50' + j^PxxxxAXi^ 



■^P(f>^(l>(l,S<P'^ + -^Px(i>(t>(i>A.Xi6(f)^ + -PxX(i>^A.Xi' 



+ ^Pxxxr^AX^'6(j). (A14) 

The next step consists in writing and solving momentum and hamiltonian contraints in 
order to integrate out the lapse and the shift functions and Ni. This was done in [15], 
where solutions for momentum and hamiltonian constraint were provided for a general 
Lagrangian. 

The constraint equations are 

V, [N~' [E] - 6]E)] = N-'Px [0 - N'dKpj 9,0, (A15) 

jl{3) + 2P - 4PxX - iV~^ {EijE'^ - E^) - 2P^xh'^di(f)dj(l) = 0. (A16) 

The action to a given order n only requires the constraint equations to be solved up to 
order n — 2. Therefore we will solve the constraints to second order in the metric and 
scalar field fluctuations. Let's employ the expansions 

a = ai + a2, (A17) 
A = + /32i, (A18) 
6 = 61 + 62. (A19) 

The momentum constraint to first order is as follows 

^did'N, - — , 
2a2 ^ 2a2 



2Hd,ai - ^did'N, - -^d,d,N' + a-^d, (dm,) = P,x<pd,6c^, (A20) 



where d' = and A^^ = S'^j. 

This can be solved to get ai: deriving both sides by and using the divergenceless 
condition for /3 the result is 

Going one step behind and using the solution found for ai, we find d^diPij. The solution 
Pij = can be chosen. From now on we will define Pi = (32i for simplicity. 
The momentum constraint to second order is 

2Hdja2 - ^d^Pj - ^Haid^ai - ^diai (^jS^^i - 6'^didj6i) (A22) 
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where 9^ = 5'^Widj. 
The solutions are 



and 



a 



2 1 



^2 



2 2Ha'^ 



P 



2H 



AH 



H 



(A23) 



(A24) 



50 
— tti 



for the vector mode. 

Let's move now to the hamiltonian constraint. To first order it is solved by 
AH 



AXP 



X 



To second order 



AH 



2PJ(\) - SP^xxX^ I ^ - tti I - 4XPx<^(50 - 12H^ai. (A25) 



2a,)[AXPx^ + 20PxxX^^ + 2XPx^6^ + SPxxxX''^ 

5'^ diOidjdct) 



^p,xx,xH<p + \d-e,] - ^^^'^ + 



X (-6/7^ + 2XPx + 4X2p^;,) [3a? - 2«2] + 4«? (SX^P^^ + 2X^P^xxx) 
+ ^ [2XPx + leX^P;,;, + SX'P^xxx] + ^ [AXP^x^ + 8XPxx0] 



AX^P^xx , 2P 



1 

+ 4 



(A26) 



We are now ready to expand the action (66) up to fourth order in perturbation theory. 
The result is provided in Appendix B. 
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Appendix B. Complete expressions for the action up to fourth order in 
perturbation theory 



The action to third order is given by 

1 
2 

where 



^(3) 



d^xdrjo^ 



-\- L^i -\- Lint 



(Bl) 



= Px 
+ P,xx 
+ ( di5(l) 



- ai5(t? + 2(t)alt 



P 



X 



ai[di5<f)J - 25<i)di5<t)diei 
P,xx 



1 



+ P,xxx 



di5(j) 



+ 



,x<t> 



3^ ^ 

6(j)(di6(l)] - 2<p6(l)di6(j)d,ei 



P,x4>aiH' + -P, 



P 



X4> 



P: 



X<j)<f> 



+ P,XX(p 



AH 



H'yu^'jf^j'yji^ 

Lint = - ^o^i'jij'jij - -^^ijdi'jjkdkdi - -^aidi-fjkdi'jjk + ^ai'jijdidjOi 

2H 1 1 2H 

T^ijljkdkdiOi + —aiiijdidjOi + -i^iijijkdkdiOi H :j-lijdajkdkOi 



(B2) 
(B3) 



+ 



dk9i 



a" 

2 „„„„„„ 2 „„„„o„ 1 



\dk9idjdk9i + -^'jijdidjOid 9i + -^di'jjkdi9idjdk9i 



a 



-^di9idj-iikdjdk9i. 



(B4) 



The complete expression for the fourth order action is 



^(4) = j £xdri 



- Haij^Y^jjk - ^I'^i'^^iikiji - H^iY^ijiki - ^liY^ijiki 

+ PxxH ^0'^' - + y + Pxh<P - y 
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Pxxxx (^^50^ - y + y - y + 

P'^^^'da.k 4:Haid^9, SHa^a^d^e^ l^l'^'ikid^ei AH ^^^^ 
—2^^ + + 2^^ ^"^^ 



a? 2a2 a2 2a^ 

5'^a2da'''d,-iki 2H^i^>^'d%dj^,i ili'^'d^d.^ki alid^Oif 
Ao? a2 2a2 



a? 

ai%d'(3^ l^i^k^^P' ifijkd^P' d,l3,d^l3' 8,13^8^13^ 
a2 2a2 ^ 2a2 + ^ 2a4 

+ + 4^^ + ^-^^'^' 



2 ^ ^ 2a2 a- 



+ 05050^ - 40'ai(5050 + 50^a?5050 - 2ifa2^c\)b^ - 20^a?50 

,a^(50<9^(50 (j)^aiS(l)Sijd'S(l)d^S(l) 



+ 20 aia2S(f) 



2(t)'^5(t)5(t)5,jd'5(i)dWi 2(P^6(f)6ijd'6(f)dWi\ f(j)'^6(j)^ 80=^ai50^ 
5 1 5 + Pxxx ' 

706q,4 

+ 50'^ai(50^ - 0^0250^ - 405a'^50 + 05aia2(50 + — — ^ - 0^aia2 

6 

<p^6<p^6ijd'6(j)d^6(f) <j)^aAjd'6(f)d^6(j) ^^al6ijd'6(j)d^6(l) 
2^2 + -2 + ^2 

'^^^^Sijd'S(j)d^6(f) 2^^aiS^Sijd'6(pdWi ^^ajSijd'Scpd^ei 



0? a2 a2 

+ Px^ - «i5050^ + 20a^5050 - 20025050 - 0^ai50 + 20^aia250 

_ 20505,^9^50/3^' _ ai505y9^509^'50 507^^-9' 509^' 50 

a2 a2 a2 

_ 26(f)6^6ijd'6(f)dWi ^ 2<Pai6(j)6,jd'6(f)dWi ^ 20507^^9^509^^1 
a2 a2 a2 

20505ij(9^509-''e2A iHaaijd'd^ei 4Ha2-fijd'dWi 



+ 



a2 y a2 a2 



a2 a2 a2 a2 

aldidjOid'dWi _ a2didjeid'd^ei _ 2ai9^9jgi9^fflgi 

4gai7,,9-9^g2 2i77f7,fc9-9^g2 aiiijd'd^92 ihjkd'd^02 
9 9 9 9 
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+ 



2a^ 

2'yij8'f3^8^6i 2-iij8' 8^618'^ 62 8ajk8'f3^8''6i _ 8jj,k8'/3W''6i 
^ ^ 

, 8^f3'8k-fij8''6i 2jjk8^f3'8''8i6i , 2aa,k8'8'6i8''8i6i 



l\lki8'8^6i8^8A A^jk8'8%8''8i62 2^ik8^ P'8''8j6^ 

^ 

I3'8ajk8''8^6, 2aajk8^6i8''8^6, ^l^ki8^6,8''8^6, 
ai8ajk8%8^8''6, -i\8aui8%8W^6^ 8ajk8%8W% 
2P'8,^,k8^8^6^ 2ai8%8j-f,k8''8Wi 28%8,-f,k8W% 



+ 



-fi8%8,-fu8'8Wi 2^. jk8^ 6,8^8^62 8a,k8%8^8^62 

28%8,-f,k8W'6, 2-f'''8%8^6,8n,k 2^i8%8W%8aik 
a4 a4 a4 

2H-fij'''8%8n,k lii'''8%8n,k ll8'6,8W%8ajk 
a2 2a2 ^ a4 

^^^'g'gig^-^igfc7,*g,7,, g-gig,-7Hg^'gig'7f g'gig^'gi.9fc7,,g'7f 
2a4 2a^ 

+ Px — 02^0^ — 20af(50 + A(paia2S(p + — 30^a^a2 

•2 2 26(f)6ij(3W^6(j) 2<i)ai5i.^'d^6(j) 207i,-/3^9^50 

+ "2 + — ^1 + + — H — 

a"^ 

a2Sij8''6(j)8^5(f) ai'yij8^5(j)8^6(t) 5^^'^ik'^ji8''54)8^5(t) 

0? a? 20? 

2ai6^6ij8'6(f)8^6(p 2S^-fij8'6(f)8Wi 2^al6ij8'6(j)8^6i 



H 1 _ _^ 

a? a? 

^ 2^a25ij8'5(t)8^6i _ 2^ai-fij8' 6^8^61 _ ^5^^^ikiji8'5(t)8^6i 

a? a? 

2S^6ij8'6(j)8W2 2<pai6ij8'6(f)8^62 2^-fij8'8W2 

a? a? a? 

, 5,k5,i8'5(l)&5ct>8^6,8'6,\ , ^ (5^" ■ , 17^^al6^^ 
+ I + Pxx [—- 30ai()0 + 

- 402^250^ - 9^^al64> + 104>^aia2S4> + ^ ^ - 6(/.^a?a2 + ^^^a 

_ <i)^6^6ij^'8^d(f) 2(j)^ai6ij(3'8^6(j) 6^^6ij8'6(f)8^6(j) 

0? 2o? 
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+ 



+ 



(i)aAjd'5<pd^<p (p5(p5ijd'5(t)d^5<p <p^al5ijd'5<pd^ 



2a2 



^'^a25ijd'505(p ^'^ai-fijd'6(pd^6(f) 3^6^6ijd'6(j)dWi 



S<i)'ai5(t)5ijd'5<f)dWi 2<p'5(hij5ijd'5(t)dWi h(p^a{5ijd'5<i)dWi 



+ 
+ 

+ 



2,p^a2Sijd'S(j)d^9i 2ai-fijd'S(l)d^9i^^ 2^'^6^6ijd'6^dW2 



,jd'5(pd^92 



4a4 



%d%\ ^ ^,k^,id^d^9,d'd^9, _ -f,j^kid'dW,d'd^9i 



(B5) 



Appendix C. Details about the computation of two- vertex diagrams 

We will provide a summary of the main steps of how to compute the two- vertex diagrams. 
We first expand the IN-IN formula as in Eq. (61) 
\2 , r rv* 



T 



T 



S^kMkA^ + B + C + D 



where 



B 
C 
D 



oo 



d^x dr] a^H^ 



int 1 



d^x a^dr] H.^ 



int 



oo 

■n* 



d^x drj a^Hj 



int ! 



d^x drj a^H,^ 



int 



oo 

■n* 



d^x dr] a^Hi 



int 1 



d^x dr] a^H- 



int 



oo 



d^x dr] a^Hl 



inf 



(CI) 

(C2) 
(C3) 
(C4) 
(C5) 



The plus and minus signs arise from having both time-ordering and anti time-ordering 
operators and give rise to different Feynman contraction rules for the field operators. 
After expanding the field operators in Hl^^ in their Fourier modes, contractions can 
be made among creation and annihilation operators. We can then integrate over space 
{d^x) using the 5 function properties and, for each diagram, we are eventually left with 
an integral over the internal momentum running in the loop combined with a double 
time-integration. Before integrating in time, it can be convenient to perform a change 
of variables x = —kr]. As to the momentum integrals, it can be helpful to consider 
the following scheme that takes into account the momentum dependence of both the 
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integrand function and the external legs as well as every momentum factor expected in 
the final expression: 

• every eigenfunction {5(f)q or 7^) generates a (2g)~^/^ factor (see Eqs. (62)- (63)); 

• spatial derivatives di and respectively produce a and a factors. No 
i factors or minus signs need to be included since the spatial derivatives always 
appear in an even number in our vertices; 

• temporal derivatives 6(f)q or 7^ originates a —{q/kY factor, e.g.: 
- e^^^m%2 (C6) 



(2g)3/2 \kJ 
where x = —kcgf]', 

• every vertex present in a given diagram provides an extra factor 

/ "'^^^ - / - 1 1^'-*^')' - - 1 T^**^-) • (C7) 

• all a" factors present in the verteces Vi contribute with extra /c" factor; 

• the polarization tensors are also responsible for producing a momentum dependence: 
they combine with one another according to the rules listes in Appendix C, Eqs. (El) 
through (E4). 

The final results for the leading two- vertex diagrams with scalar loops are 
Gaaix^Faaik) = ^ (-4 - Ax*^ - 25x*^ + x*^) lu (A/H,) , (C8) 

Gabix*)Fatik) = (19 - 71x*^ + 41x*^) In (A/if,) , (C9) 

Gtb{x*)F,,{k) = (32 + 384x*2 + 87x*^ + 4x*6) In (A/H,) , (CIO) 

where the functions G and F were introduced in Eq. (92). 



Appendix D. Details about the computation of one-vertex diagrams 

Let us now consider in details the contributions and analytic results from one-vertex 
diagrams. We will first recall the remark made at the end of Sec. 4 about the one- 
vertex diagrams with tensor loops. This is important in order to quickly recognize the 
diagrams that provide contributions that can be absorbed in a renormalization constant 
and whose computation can be therefore avoided. Indeed, the loop contributions that 
are interesting are those that originate from a momentum loop integral with a real 
dependence from the external momentum, i.e. that result in a logarithmic function of 
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momentum in the final result. On the other hand, a four- leg vertex that is responsible 
for a diagram with the following mathematical form 

r+oo 

Pol{x*) / (Dl) 



(where Pol{x*) is a polynomial function of the horizon crossing time x* = —kt]*, a is 
a positive integer), produces a disconnected-like contribution to the two-point function 
which can be absorbed into a left-over renormalization constant. For this reason, tensor 
vertices have not been included in Eq. (71) and have not been taken into account in 
our one-vertex diagrams computations: by looking at the complete expression for the 
Lagrangian provided in Appendix B (Eq. (B5)), it is easy to verify that all of the 
graviton non zero contributions have a form as in Eq. (Dl). For the same reason, the 
only leading order scalar diagrams that we have taken into account as far as the one-loop 
contributions are concerned, are the ones from the vertices given in Eqs. (108) through 
(112), which we rewrite below 

V(o) = "2^0 , (D2) 

= 24|^ (aj,) , (D3) 

AS 

1^(2) = -12— (D4) 

^(3) = TTr^2 {Sijdi6(t>dj6(f)) , (D5) 

V^(4) = -6^^505,, {dM) (djO^) . (D6) 
These can be rewritten as 

^(1^) ~ didjd"^ (5050) didjd"^ (d'^6(f)6^^ , (D7) 

V(iB) ~ didjd~^ (5050) d^djd-^ (^5imdiS(j)dm5<i)) , (D8) 

V^,c) - didjd~^ (5050) (9,9"2 (5,<5050) , (D9) 

V(3A) ^ d-^ (9^5050) SijdMdj6(j), (DIO) 

V(3B) ^ d-^ (^6ijdi64>dj6^^ 6imdi64>dmS(p, (Dll) 

V(4) ~ 6^6ijdi6(f)djd-^ (5050) . (D12) 

As to V(o) and V(2), using the definition of A and S from Eqs. (11) and (12) and referring 
back to Table 1, it is easy to verify that they cancel each other to leading order in the 
slow- variation parameters. 

From the expansion of the Schwinger-Keldysh formula (61) we have 

50,-;). D z 1 ( [ifL1(7^'), 50,-;50,-;] ) . (D13) 



35 



Before proceeding with the derivation of the analytic expressions of the diagrams 
corresponding to V(i), V(3) and V(4), it is important to notice that, for each one of 
them, both "regular" and "singular" contractions exists. To understand what we mean, 
let us consider the diagrams originating from V(i^) ~ didjd~^ ^5050^ didjd~'^ (^d^ScpScj)^ . 
If we call qi, g*2, ^3 and respectively the momenta associated with the four internal 
legs, we count three possible different types of contractions, i.e. 

kiq2 X hqi X 93^4, (D14) 

hqi X k2q3 X ga^i, (D15) 

kiq4 X k2qi x ^2^3, (D16) 

where the " symbol stands for a contraction between the two wavefunctions to which the 
momenta are associated. It is easy to check that the second and the third permutations 
respectively provide a factor 

Q-ki] ■ (q + h]] 

(D17) 



^2 



\q- ki\'^\q + k2\'^ 



q + ki 



\q- k2\'^\q + ki\ 



(D18) 



inside the momentum integral, where q is the momentum running inside the loop. The 
first permutation instead gives a factor 

1 2 



q^ 



ki + k2) ■ {q- q) 



\ki + k2\^^ 



(D19) 



which is apparently singular, also considering momentum conservation. 

We will now show how these "singular" diagrams end up not contributing to the final 

result at all. 



From Via we get a singular contribution 

2 



d^q 



2 / dxTJx)-^ 



{ki + k2) -{q-q} 



4^4 



k q 



\ki + /C2I \q- q\' 



2 / dxTJx) 



{ki + k2) ■ (g - q) k'^q^ 



— * — * t r 

\ki + /C2I \q- q\' 



(D20) 



(where we have condensed the products of wavefunctions inside the symbols and T},) 
which is equal and opposite to the one from Vib- Similarly, we verified that the singular 
diagrams from V^a 



d^q 



2 / dxTJx] 



\ki + k2\ 



2 / dxTAx 



q\ki ■ k2) 

\<i-q\ 



(D21) 
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(11,(11 




Figure 5. Diagrammatic representation of the one-vertex interaction diagrams. 
The dots on each end of a hne indicate temporal derivatives of the corresponding 
wavefunctions. 



cancels the ones from V^b- Finally we can show that the singular contributions from 
Vic and V4 (surviving after we sum up over the permutations) 

{h + h) ■ (<? - q)] \q- ih + h)] kS^' 



ic 



Va 



^ ^ 2 2 
1^1 + k2\ \q- q\ 



k'^q^q- {ki + /C2) 



(D22) 
(D23) 

(D24) 



where we define e = ki + k2- From Eq.(D24) it is apparent that the integral goes to zero 
for symmetry reasons. The same argument applies to Vic- Let us write it as 

^■^^^^-^^ (D25) 



\ki + k2 

give a zero contribution as well. As to V4, one can write it as 
[ r 

a qq q 



d'^qq 



where 6 = q — q and we assumed that 5 is independent of the integration variable q and 
so can be taken outside the integral. 

We are now left with dealing with the "regular" diagrams. The time integrals are 
diagrammatically represented in Fig. (5). The momentum integrals can be easily 
performed. We report the final results 



Gi{x*)Fi{k) = ^ (-5 - 5a:*2 - ISx*^) i^^^^, 
G3(x*)F3(A;) = ^ (-405 - 745X*' - 366x*^) ^"^Z^* 



fc3 



(D26) 
(D27) 



37 



G^{x*)F^{k) = ^ (5 + 5x*2 + 18x*'^) h^^^/lA, (D28) 
The ln(A/i7*) terms come from ultraviolet divergences. 

Appendix E. Polarization tensor equations 

The following relations can be useful in the calculations involving the tensor modes 

k"^ sin^ e 

k'k^eiiq)sl^iq) = 2esm'e, (E2) 
{k^-{q-z?? 



e^^mk^ = (El) 



{k'-{q-zyf{-k^ + {q+-'' 



(E4) 

— * 

where z = \q — k\ and sums are taken over repeated (spatial and polarization) indices. 
These equations can be easily derived as follows. We can choose a spatial coordinate 
frame so that the vector k points along the third spatial direction, then the components 
of eij{k) are en = —622 = l/\/2, £12 = £21 = ^i/V^, = £31 = (the plus and minus 
signs refer to the two polarization states of the tensor). The tensor eij{q) is obtained 
by rotation along the direction of g = (sin ^ cos 0, sin ^ sin 0, cos ^). The result is 

cos^ 9 1 sin^ 9 

= '" = "71' = 

i cos 9 — sin 9 cos 9 



£12 — £21 — ± — 7=^, £^13 — £^31 



£23 = ^32 = ^^yf- (^^) 

where we set = 0. 

Similarly, eij{z) is obtained by rotating eij{z) along the direction of z = (sin ^2, 0, cos ^2), 
where sin6'2 = sin^(g/^) and cos^^ = {qcos9 — k)/z. 

Using the matrices ejj(g) and eij{z) thus constructed and considering the polarization 
tensor orthogonality condition q^£ij{q) = 0, Eq.(El) through (E4) can be 
straightforwadly derived. 
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